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Abstract 

The critical surface for random-cluster model with cluster-weight g > 4 on iso¬ 
radial graphs is identihed using parafermionic observables. Correlations are also 
shown to decay exponentially fast in the subcritical regime. While this result is 
restricted to random-cluster models with q > A, \i extends the recent theorem of 
[6] to a large class of planar graphs. In particular, the anisotropic random-cluster 
model on the square lattice is shown to be critical if = g, where and 

Ph denote the horizontal and vertical edge-weights respectively. We also mention 
consequences for Potts models. 


1 Introduction 

1.1 Motivation 

Random-cluster models are dependent percolation models introduced by Fortuin and 
Kasteleyn in 1969 [25]. They have been an important tool in the study of phase transi¬ 
tions. Among other applications, the spin correlations of Potts models get rephrased as 
cluster connectivity properties of their random-cluster representations, which allows for 
the use of geometric techniques, thus leading to several important applications. Never¬ 
theless, only few aspects of the random-cluster models are known in full generality. 

The random-cluster model on a finite connected graph G = (V[G'],T[G]) is a model 
on edges of this graph, each one being either closed or open. A configuration can be seen 
as a random graph, whose vertex set is V[G] and whose edge set is the collection of all 
open edges; a cluster is a connected component for this random graph. The probability 
of a configuration is proportional to 

n n (i-p.)'T 

e open e closed 

where the edge-weights pe G [0, 1] (for every e G £[G]) and the cluster-weight q G (0, oo) 
are the parameters of the model. Extensive literature exists concerning these models; we 
refer the interested reader to the monograph of Grimmett [27] and references therein. 

For q > 1, the model can be extended to infinite-volume lattices where it exhibits a 
phase transition. In general, there are no conjectures for the value of the critical surface, 
i.e. the set of (pe)ee£[G] for which the model is critical. In the case of planar graphs, there 
is a connection (related to the Kramers-Wannier duality [37, 38] for the Ising model) 
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between random-cluster models on a graph and on its dual with the same cluster-weight 
q and appropriately related edge-weights. This relation leads, in the particular case of 
7? (which is isomorphic to its dual) with pe = p for every e, to the following natural 
prediction: the critical point Pc{q) is the same as the so-called self-dual point, which has 
a known value y/q /(1 -|- -^/q). This was proved recently in [6] for any q > 1 (see also [20]). 
Furthermore, the size of the cluster at the origin was proved to have exponential decay 
tail if p < Pc{q)- 

The critical point was previously known in three famous cases. For q = 1, the model 
is simply Bernoulli bond-percolation, proved by Kesten [35] to be critical at Pc(l) = 1/2. 
For q = 2, the self-dual value corresponds to the critical temperature of the Ising model, 
as first derived by Onsager [43]; one can actually couple realizations of the Ising and 
random-cluster models to relate the critical points of each, see [27] and references therein 
for details. Finally, for q > 25.72, a proof is known based on the fact that the random- 
cluster model exhibits a first order phase transition; see [39, 40]. 

A general challenge in statistical physics is to understand the universal behavior, i.e. 
the behavior of a certain model, for instance the planar random-cluster model, on different 
graphs. A large class of graphs, which appeared to be central in different domains of 
planar statistical physics, is the class of isoradial graphs. An isoradial graph is a planar 
graph admitting an embedding in the plane in such a way that every face is inscribed in 
a circle of radius one. In such a case, we will say that the embedding itself is isoradial. 



Figure 1: The black graph is the isoradial graph. White vertices are the 
vertices of the dual graph. All faces can be put into a circumcircle of radius 
one. Dual vertices have been drawn in such a way that they are the centers of 
these circles. 

Isoradial graphs were introduced by Dufhn in [17] in the context of discrete complex 
analysis. The author noticed that isoradial embeddings form a large class of embeddings 
for which a discrete notion of Cauchy-Riemann equations is available. Isoradial graphs 
first appeared in the physics literature in the work of Baxter [1], where they are called 
Z-invariant graphs; the so-called star-triangle transformation was then used to relate 
the free energy of the eight-vertex and Ising models between different such graphs. In 
Baxter’s work, Z-invariant graphs are obtained as intersections of lines in the plane, and 
are not embedded in the isoradial way. The term isoradial was only coined later by 
Kenyon, who studied discrete complex analysis and the graph structure of these graphs 
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[34], Since then, isoradial graphs were used extensively, and we refer to [9, 14, 30, 34, 42] 
for literature on the subject. 

In the present article, we study the random-cluster model on isoradial graphs. 

1.2 Statement of the results 

All the graphs which we will consider in this paper will be assumed to be periodic^ in the 
sense that they will carry an action of the square lattice 1? with finitely many orbits; 
indeed, this is often a crucial hypothesis in the usual arguments of statistical mechanics. 
Nevertheless, some of our results extend to a more general family of isoradial graphs, 
which is why we first introduce the following, weaker condition. 

Let e be an edge of an isoradial embedding: it subtends an angle 0^ G (0, vr) at the 
center of the circle corresponding to any of the two faces bordered by e; see Fig. 1. Fix 
0 > 0, and let Goo = (V[Goo], ^^[^* 00 ]) be an infinite isoradial graph. The graph is said to 
satisfy the bounded-angle property if the following condition holds: 

(BAPe) For any e G 8[Goo\^ 9 < Og < n — 9. 

In order to study the phase transition, we parametrize random-cluster measures with 
cluster-weight g > 4 with the help of an additional parameter fd > 0. For (3 > 0, define 
the edge-weight Pei/3) G [0,1] for e G ^[ 6 * 00 ] by the formula 

Peil3) _ sinh[^^^^^] 

[^ - Pei/3)] y/q sinh[^] 

where the spin a is given by the relation 



The infinite-volume measure on Goo with cluster-weight g > 4, edge-weights (pe(/3) : e G 
T[Goo]) and free boundary conditions (see next section for a formal definition) is denoted 
by /3,9- 

Remark 1 In the case of the square lattice, one gets Pe{/3) = This does not 

quite match what one obtains in the setup of the coupling between the Potts and random- 
cluster models: the bond-parameter corresponding to the q-state Potts model at inverse 
temperature /d is equal to 1 — This simply means that what we will denote here by 

[3 should not be interpreted as an inverse temperature as such, but simply as a parameter 
according to which a phase transition can be defined. 

Let I ■ I be the Euclidean norm. 

Theorem 2 Let q > A, 9 > 0 and /3 < 1. There exists c = cif3,q,9) > 0 such that for 
any infinite isoradial graph Goo satisfying (BAP 0 ), 

is connected to v by an open path) < exp[—c|a; — y\], 

for any u, n G Goo • 

This theorem implies that the edge-weights Pe = Pe{3) are critical in the following 
sense. 
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Theorem 3 Let q>A, 0 > 0. For any periodic isoradial graph Goo-' 

1. The infinite-volume measure is unique whenever [3^1. 

2. For (3 < 1, there is 4 >g^ ^ almost surely no infinite-cluster. 

3. For f3 > 1, there is 4 >Goo yq-^^f^ost surely a unique infinite-cluster. 

In fact, what we will prove is the following, slightly weaker result in the more general 
setup of graph satisfying the bounded-angle property: 

Theorem 4 Let q > 4, 6 > 0. For any infinite isoradial graph Goo satisfying (BAPe): 

1. For f3 <1, there is fcoo yq-^^^T^ost surely no infinite-cluster. 

2. For [3 > 1, there is ^ q-almost surely a unique infinite-cluster. 

It will be shown that in the periodic case, or in any case for which the set 'Dq^Coa 
of f3 such that there are more than one infinite-volume random-cluster measure is of 
everywhere dense complement (see Proposition 8 below). Theorem 4 implies Theorem 3. 
Since this will be the only place where periodicity will be used, most statements of this 
article are phrased (and proved) in the more general bounded-angle setup. 

The theorems were previously known for two specific choices of q: when g = 2, 
the model was proved to be conformally invariant when /3 = 1 in [15], thus implying 
the different theorems; for percolation (i.e. the case q = 1), Manolescu and Grimmett 
[28, 29, 30] showed the corresponding statements (these papers contain more delicate 
properties of the critical phase /? = 1 as well). 

The main tools used in the proofs are the so-called parafermionic observables. These 
observables were first introduced in [44] for critical random-cluster models on If with 
parameter q G [0,4], as (anti)-holomorphic parafermions of fractional spin a G [0,1], given 
by certain vertex operators. So far, holomorphicity was rigorously proved only for q = 2 
(which corresponds to the Ising model) and probably holds exactly only for this value. 
In this case, the observable can be used to understand many properties on the model, 
including conformal invariance of the observable [15, 44] and loops [11, 33], correlations 
[12, 13, 32] and crossing probabilities [8, 10, 19]. Inspired by similar considerations, one 
can also compute the critical surface of any bi-periodic graph [41, 16]. 

Our proof uses an appropriate generalization of these vertex operators to random- 
cluster models with g > 4. Even though exact holomorphicity is not available, the 
observable can still be used efficiently. Interestingly, the spin variable becomes purely 
imaginary and does not possess an immediate physical interpretation. However, this 
allows us to write better estimates even in the absence of exact holomorphicity. It also 
simplihes the relation between our observables and the connectivity properties of the 
model. 

For /3 7 ^ 1, we prove that the observables behave like massive harmonic functions and 
decay exponentially fast with respect to the distance to the boundary of the domain. 
Translated into connectivity properties, this implies the sharpness of the phase transition 
at /d = 1. 

The fact that isoradial graphs are perfect candidate for constructing parafermionic 
observables is reminiscent from both the works of Duffin and Baxter. Indeed, these 
works highlighted the fact that isoradial graphs constitute a general class of graph on 
which discrete complex analysis and statistical physics can be studied precisely. 
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Application to inhomogeneous models The inhomogeneous random-cluster models 
on the square, the triangular and the hexagonal lattices can be seen as random-cluster 
models on periodic isoradial graphs. Theorem 3 therefore implies the following: 

Corollary 5 The inhomogeneous random-eluster model with cluster-weight q > A on the 
square, triangular and hexagonal lattices 1?, T and El have the following critical surfaces: 

Pi P2 
1i 

1 - Pi 1 - P2 

Pi P2 P3 , Pi P2 , Pi P3 , P2 P3 
i-i-i-^ i-i-^ 1 -i-^ l -i- “ 

1 - Pi 1 - P2 1 - P3 1 - Pi 1 - P2 1 - Pi 1 - P3 1 - P2 1 - P3 
Pi P2 P3 _ Pi , P2 , P3 , 2 

1 - Pi 1 - P2 1 - P3 1 - Pi 1 - P2 1 - P3 

where pi,p 2 (resp. Pi,P 2 ,P 3 ) cire the edge-weights of the different types of edges. 

For percolation, Corollary 5 was predicted in [46] and proved in [36, Section 3.4] for the 
case of the square lattice and [26, Section 11.9] for the case of triangular and hexagonal 
lattices. 

Let us also mention that the critical parameter of the continuum random-cluster 
model can be computed using the fact that it is the limit of inhomogeneous random- 
cluster models on the square lattice with (pi,P 2 ) —)■ (0,1). We refer to [31] for a precise 
definition of the models, which are connected to the one dimensional Quantum Potts 
model. The parameters of the models are usually referred to as A, 5 > 0, where A and 
6 are the intensities of the Poisson Point Process of bridges and deaths respectively. In 
such case. Theorem 3 implies that the critical point is given by A/5 = g for g > 4. 

Application to Potts models Potts models on G with g colors and coupling constants 
(Je : e G T[G]) can be coupled with random-cluster model with cluster-weight g and edge- 
weights Pe = 1 — exp[—Jg]. As a consequence. Theorem 3 shows the following: 

Corollary 6 Let g > 4 and 6 > 0. For any infinite periodic isoradial graph Goo, the q- 
state Potts models on isoradial graphs with coupling constants — log[l — Pe(l)], e G £[Goo] 
is critical. 

1.3 Open questions 

Exact computations can be performed for the random-cluster model at criticality (see 
[2]), and despite the fact that they do not lead to fully rigorous mathematical proofs, 
they do provide insight and further conjectures on the behavior of these models at and 
near criticality. Let us mention a few open questions. 

1. Parafermionic observables were used when 1 < g < 4 to prove that the phase tran¬ 
sition is continuous [18, 21]. Moreover, it is conjectured that among all random-cluster 
models defined on planar lattices, the phase transition is of first order if and only if g 
is greater than 4. Interestingly, the parafermionic observable exhibits a very different 
behavior for g < 4 and g > 4, which raises the following question. 

Question 1. Can the change of behavior of the observable be related to the change 
of critical behavior of the random-cluster model? 


on l} 
on T 
on El 
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2. In the work [6], the critical value for the random-cluster model on the isotropic square 
lattice has been computed for any q>l. Parafermionic observables on isoradial graphs 
also make sense for g < 1 (see [18, 44]), which leads to the following question. 

Question 2. Use the parafermionic observable to compute the critical point on 
isoradial graphs (or simply on 1?) for any q G (0,4)? 

3. More generally, parafermionic observables have been found in a number of critical 
planar statistical models, see [22, 45] and references therein. They have sometimes been 
used to derive information on the models (see [18] for random-cluster models and [3, 4, 5, 
23, 24] for (9(n)-models and self-avoiding walks). There are many potential applications of 
these observables which deserve a closer look and we refer to the literature for additional 
information on these questions. 

4. As mentioned earlier, the fact that random-cluster models on 7? undergo a hrst order 
phase transition is currently known for q > 25.72; see [39, 40]. The main ingredient is the 
Pirogov-Sinai theory, which shows that the 0^2 1 ^-probability that the origin is connected 
to distance n decays exponentially fast in n. Interestingly, Grimmett and Manolescu [30] 
used the star-triangle transformation to relate probabilities of being connected to distance 
n for percolation on different isoradial graphs. From [1], the star-triangle transformation 
is known to extend to critical random-cluster models and it seems plausible that the tech¬ 
niques in [30] can be combined to results in [39, 40] to prove that the 0^^ ^-probability 
random-cluster models that the origin is connected to distance n decays exponentially fast 
in n whenever q > 25.72. This would show some kind of universal behavior: hrst order 
phase transition is common to any random-cluster model with large enough cluster-weight 
on isoradial graphs. Note that Pirogov-Sinai theory extends partially to this context (al¬ 
though likely with different bounds due to the fact that the graphs involved would have 
different combinatorics). 

Question 3. Show that random-cluster models on any isoradial graph undergo a hrst 
order phase transition when q is large enough. 

5. Let us conclude with a pair of more technical questions: How to release the periodicity 
assumption. For instance, how to show Proposition 8 for isoradial graphs satisfying only 
the bounded-angle property? Can the results be extended to (non-periodic) isoradial 
graphs which do not satisfy the bounded-angle property? 

Organization of the paper. Section 2 gives an overview of probabilistic properties of 
the random-cluster model. It also introduces the observable. Section 3 contains a deriva¬ 
tion of a representation formula, similar to the formula for massive harmonic functions, 
which is then used to provide bounds on the observable. Section 4 then contain the proof 
of Theorem 2 and Section 5 the proofs of Theorem 3 and its corollaries. 

2 Basic features of the model 

We start with an introduction to the basic features of random-cluster models. Details 
and proofs can be found in Grimmett’s book [27]. 
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Isoradial graphs As mentioned earlier, an isoradial graph G = (V[G],£^[G]) is a planar 
graph admitting an embedding in the plane in such a way that every face is inscribed in 
a circle of radius one. In such case, we will say that the embedding is isoradial. For the 
isoradial embedding, we construct the dual graph G* = (V[G*],as follows: V[G*] 
is composed of all the centers of circumcircles of faces of G. By construction, every face 
of G is associated to a dual vertex. Then, T[G*] is the set of edges between dual vertices 
corresponding to adjacent faces. Edges of £[G*] are in one-to-one correspondence with 
edges of £[G]. We denote the dual edge associated to e by e*. 

From now on, we work only on an infinite isoradial graph Goo embedded in the isoradial 
way. Note that the graph is not a priori periodic. 

Definition of the random-cluster model. The random-cluster measure can be de- 
hned on any graph. However, we will restrict ourselves in this article to the graph Goo 
and its connected finite subgraphs. Let G = (V[G], T[G]) be such a subgraph. We denote 
by dG the vertex-boundary of G, i.e. the set of sites of G linked by an edge to a site of 

Goo\G. 

A configuration ca on G is a random subgraph of G having vertex set V[G] and edge 
set included in £[G]. We will call the edges belonging to u open, the others closed. Two 
sites u and v are said to be connected, if there is an open path — i.e. a path composed 
of open edges only — connecting them. The previous event is denoted by u < —)■ v (we 
extend the notation U < — > V to the event that there exists an open path from a vertex 
of the set U to a vertex of the set V). The connected components of u will be called 
clusters. 

A set ^ of boundary conditions is given by a partition of dG. The graph obtained 
from the configuration u by identifying (or wiring) the vertices in dG that belong to the 
same component of ^ is denoted by a; U Boundary conditions should be understood as 
encoding how sites are connected outside of G. Let k{u,^) be the number of connected 
components of ca U The probability measure ficpq random-cluster model on G 

with parameters p = (pe : e G T[G]) G [0, q G (0, cx)) and boundary conditions ^ is 

defined by 

4,p.,({‘p}) = -. (2-1) 

^G,p,q 

for any subgraph u of G, where Zq p g is a normalizing constant referred to as the partition 
function. When there is no possible confusion, we will drop the reference to parameters 
in the notation. 

Three specific boundary conditions Three boundary conditions will play a special 
role in our study; 

1. Free boundary conditions are the boundary conditions obtained by the absence of 
wiring between boundary vertices. The corresponding measure is denoted by 0^ p g. 

2. Wired boundary conditions are the boundary conditions obtained by wiring every 
boundary vertices. The corresponding measure is denoted by ficpq- 

3. Assume that dG is a self-avoiding polygon in Goo, and let a and b be two sites 
of dG. The triple (G, a, b) is called a Dobrushin domain. Orienting its boundary 
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counterclockwise defines two oriented boundary arcs dab and dba', the Dobrushin 
boundary conditions are defined to be free on dab (there is no wiring between these 
sites) and wired on dba (all the boundary sites are wired together). We will refer 
to those arcs as the free and the wired arcs, respectively. The measure associated 
to these boundary conditions will be denoted by We will often use the dual 

arc d*fj adjacent to dab instead of dab- See Fig. 2. 

Remark 7 The term “Dobrushin boundary condition” usually refers to mixed +/— bound¬ 
ary condition in the setup of the king model; however the main idea is the same here, 
this choice of boundary condition forces the existence of a macroscopic boundary between 
two regions in the domain (+!— for the king model, open and dual-open in the case of 
the random-cluster model), which is why we use the same term here. 

The domain Markov property One can encode, using appropriate boundary condi¬ 
tions the influence of the conhguration outside F on the measure within it. In other 
words, given the state of edges outside a graph, the conditional measure inside F is a 
random-cluster measure with boundary conditions given by the wiring outside F. More 
formally, let G be a graph and fix F C T[G]. Let X be a random variable measurable in 
terms of edges in F (call X£:[g]\f the a-algebra generated by edges of £[G] \ F). Then, 

where f denotes boundary conditions on G, is a, conhguration outside F and is the 
wiring inherited from ^ and the edges in ip. We refer to [27, Lemma (4.13)] for details. 

Comparison between boundary conditions Random-cluster models with param¬ 
eter g > 1 are positively correlated] see [27, Theorem (2.1)]. It implies that for any 
boundary conditions 'ip < f, (meaning that the wirings existing in pj exist in p as well), 
we have 

4,p.,(- 4) < 4.py-4) (2.2) 

for any increasing event A. We immediately obtain that 0^ p ^(A) < (pQ^ g{A) < (pQ p^^{A) 
for any increasing event A and any boundary conditions f. 


Planar duality In two dimensions, one can associate to any random-cluster measure 
with parameters p and g on G a dual measure. Let us focus on the case of free and wired 
boundary conditions. 

Consider a conhguration oj sampled according to (pG,p,q- Construct an edge model on 
G* by declaring any edge of the dual graph to be open (resp. closed) if the corresponding 
edge of the primal graph is closed (resp. open) for the initial random-cluster model. The 
new model on the dual graph is then a random-cluster measure with wired boundary 
conditions and parameters p* = p*(p, g) G [0,and g satisfying 


Pe* = 


(1 -Pe)g 
(1 -Pe)q+Pe 


or equivalently 


Pe*Pe 


(1 -P*.)(l -Pe) 


= q- 


This relation is known as the planar duality. Similarly, the dual boundary conditions of 
wired boundary conditions are free boundary conditions. See [27, Section 6.1]. 
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Infinite-volume measures A probability measure 0 on is called an infinite- 

volume random-cluster measure on Gao with parameters p and q if for every event A E 
and any finite G C Gao, 

for 0-almost every ^ G where -7 ^£-[Goo\g] is the a-algebra generated by edges in Gao \ G. 

The domain Markov property and the comparison between boundary conditions allow 
us to define an infinite-volume measure as the (increasing) limit of a sequence of random- 
cluster measures in finite nested graphs Gn Gao with free boundary conditions. In such 
cases, the sequence of measures is increasing. We denote the corresponding limit mea¬ 
sure fcoopg- Similarly, one can construct the measure fcoopg considering measures 
on nested boxes with wired boundary conditions. The comparison between boundary 
conditions ( 2 . 2 ) implies that for any increasing event A 

<t’Gv.,W ^ (2.3) 

Section 4 of [27] presents a comprehensive study of this question. 

The diamond graph of a Dobrushin domain Let Gao be an infinite isoradial graph. 
Define = (V[G^], T[G^]) to be the graph with vertex set V[Goo] U V[G^] and edge 
set given by edges between a site x of V[Goo] and a dual site v of V[G^] if x belongs 
to the face corresponding to v. It is then a rhombic graph, i.e. a graph whose faces are 
rhombi; see Fig. 2. To emphasize the distinction with edges of Gao, and G'^, we refer 
to the latter as diamond edges. 

We now define the diamond graph in the case of Dobrushin domains. Let {G, a, b) be 
a Dobrushin domain. The diamond graph G^ = (VfG*],£^[G^]) is the subgraph of G^ 
composed of sites in V[G] U V[G*] U and of diamond edges between them; see Fig. 2 
again. 

Loop representation on a Dobrushin domain Let [G, a, b) be a Dobrushin domain. 
In this paragraph, we aim for the construction of the loop representation of the random- 
cluster model. 

Consider a configuration u, it defines clusters in G and dual clusters in G*. Through 
every face of the diamond graph passes either an open edge of G or a dual open edge of 
G*. Therefore, there is a unique way to draw Eulerian {i.e. using every edge exactly once) 
loops on the diamond graph — interfaces, separating clusters from dual clusters. Namely, 
loops pass through the center of diamond edges, and in a face of the diamond graph, loops 
always makes a turn so as not to cross the open or dual open edge through this face; see 
Figure 2. We further require that loops cross each diamond edge orthogonally. Besides 
loops, the configuration will have a single curve joining the edges adjacent to a and b, 
which are the diamond edges Ca and connecting a site of dab to a dual site of This 
curve is called the exploration path] we will denote it by 7 . It corresponds to the interface 
between the cluster connected to the wired arc dba and the dual cluster connected to the 
free arc d*^^. 

This provides us with a bijection between random-cluster configurations on G and 
Eulerian loop configurations on G'^. This bijection is called the loop representation of the 
random-cluster model. We orientate loops in such a way that they cross every diamond 
edge e with the end-point of e in V[G] on its left, and the end-point of e in V[G*] on its 
right. 
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Let p G (0, The probability measure can be nicely rewritten (using Euler’s 

formula) in terms of the loop picture: 



({^}) 



eGo; 


where is a normalizing constant and x = [xe : e G T[G]) G (0, is given by 

Pe 

(l-Pe)v^- 



Figure 2: Construction of the diamond graph and the loop representation. 


Critical weights for isoradial graphs In the case of isoradial graphs, a natural family 
of weights can be defined. Let 

_ Pe(l) _ sin[^] 

(1-Pe(l))v^ sin[^^^]' 

The bounded angle property immediately implies that weights are uniformly bounded 
away from 0 and 1. 

This family of weights on isoradial graph is self-dual, in the sense that the dual of a 
random-cluster model with edge-weights x = (xg : e G T[G]) is a random cluster model 
on the (isoradial) dual graph with edge-weights (xg* : e* G £[G*]). 

Let us stress out the fact that many other families of weights x are self-dual. Never¬ 
theless, this family will play a special role for reasons that will become apparent later in 
the article. 

Fix /d > 0. From now on, we will consider random-cluster measures 

■^G,/3x,ij eSo; 

for any conhguration u. Note that (0G,/3,g)* = 4 * /3-i q- 
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Phase transition and critical point in the periodic case In this paragraph, iso- 
radial graphs are assumed to be periodic. We aim to study the behavior of (j^cpq ''^hen 
/3 varies from 0 to oo. Positive association of the model implies that 

for any increasing event A and (3 < (3' (in such case, (i)% ^ q{A) is said to be stochastically 
dominated by 0^^, ^(A)). The previous inequality extends to the inhnite volume. It is 
therefore possible to dehne 

Pc = sup{^ > 0 : 0Goc,/3,g(O ^^ oo) = 

where 0 <—> oo denotes the fact that 0 is contained in an inhnite open path. This value 
is called the critical point. 

We used I3 g dehne the critical point but the inhnite-volume measure is not 
necessarily unique. Nevertheless, it can be shown that for a hxed q > 1 and a periodic 
isoradial graph Goo, uniqueness can fail only on a countable set 'Dq c^- More precisely: 

Proposition 8 Let G^o be a periodic isoradial graph. There exists an at most countable 
set Vq^Goo C (0, oo) such that for any P p. T)q^G^, there exists a unique infinite-volume 
measure on Goo with parameters P and q>l. 


Proof The proof follows the argument of [27, Theorem (4.63)] quite closely, so we only 
give a sketch here. Dehne the free energy per unit volume in a hnite box as 

We have 

d 1 

^^g(/^) = Je\^\ Yl ^ [0,max{a:e : e G ^[Goo]}]; 


in particular, Hq is convex. Now, let G increase to cover the whole lattice. A classical 
argument of boundary-area energy comparison (see the lines following [27, Equation 
(4.71)1) shows that the limit H{P) of Hq[P) exists and does not depend on the boundary 
condition f. 

Since 77 is a uniform limit of convex functions, it is convex, and therefore diherentiable 
outside an at most countable set 'Dq c^- We may check that for any P p. both 

^Hq(P) and -^Hq{P) converge to the same limit, which is also equal to -^H{P). Hence, 
for any such /d. 


lim , Xed'r’fl „(e is open) 

' ^ e&E[G] 


lim , Xgd*)-« „(e is open) 

G^G^ EG ^ ^ > 

' ^ e&E[G] 


which is enough to guarantee that the measures 0 Goo,/?,<? ^G^,i 3 ,q coincide; this in turn 

implies uniqueness of the Gibbs measure for all P ^ T)q^G^- D 


Since the inhnite-volume measure is unique for almost every P (at hxed q), for any 
inhnite-volume measure pGco,i3,qj 


= 0 A P < Pc 
> 0 A P > Pc 


0GooA 9(O ^^ oo) 
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(2.4) 









Observables for Dobrushin domains Fix a Dobrushin domain {G, a, b) and consider 
the loop representation of the random-cluster model. Following [44], we now define an 
observable F on the edges of its diamond graph, i.e. a function F : £[G^] —)■ M+. Roughly 
speaking, F is a modification of the probability that the exploration path passes through 
the center of an edge. First, we introduce the following definition; the winding Wr(;s, z') 
of a curve F between two edges 2 ; and z' of the diamond graph is the total rotation (in 
radians) that the curve makes from the center of the diamond edge 2 ; to the center of the 
diamond edge z'. 

Let g > 4. We define the observable F for any diamond edge e G £[G^] by 
where 7 is the exploration path and a > 0 is given by the relation 



For cr > 0 to exist, q needs to be larger than 4, hence the hypothesis in the theorem. We 
define the function F by 

ne) = (e-'W-l'''**!.,,) . (2J) 

Remark 9 The observable G{e) = 12^,07 , where sin((T7r/2) = ^/q/2, was 

introduced in the case q < A in [44] for the square lattice. When weights are critical, one 
obtains around each vertex v 

G{NW) - G{SE) = i[G(iVF) - G{SW)], 

where NW, SE, NE and SW are the your edges incident to v indexed in the obvious 
way. This relation can he seen as a discrete version of the Cauchy-Riemann equation. The 
observable is then a holomorphic parafermion of spin a (which itself is a real number in 
[0, V\). For q > 4, a is purely imaginary and does not have an obvious physical meaning; 
it would nonetheless he interesting to find one. In this article, one could work with the 
corresponding G for q > A, but the definitions in (2.5) and (2.7) are easier to handle for 
the applications we have in mind. 


3 A representation formula for the observable 


Let (G, a, b) be a Dobrushin domain. In this section, we estimate the sum of F over a 
set E C £[G^] in various ways. Let be the set of inner faces of G^. Any / G is 
bordered by four edges in £[G^], which we label counterclockwise A, B, G and D, so that 
locally the loops (or the exploration path) go from / to the outside when crossing A and 
G, and from the outside to / when crossing B and D] see Figure 3. There are a priori 
two ways to do so, but the choice will be irrelevant. 


Lemma 10 Fix /? > 0 and q > A. For every face f G 


F{B) + F{D) = A,(Px,) |F(2l) + F{C )], 


(3.1) 


where e is the edge ofG passing through f, and Ag is given by Ae(x) = e 




X -f e ""a 


A similar statement was used in [7] to derive massive harmonicity of the observable 
when g = 2 on the square lattice. This enabled to compute the correlation length of the 
high temperature Ising model. Observe that A(xe) = 1. 
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Proof Consider the involution s on the space of conhgurations which switches the state 
(open or closed) of the edge of G passing through /. 

Let e be an edge of the diamond graph and denote by the 

contribution of oj to F{e) (here p{uj) is the probability of the conhguration u). Since s is 
an involution, the following relation holds: 

F{g) — ^ ^ Coj — 2 ^ ^ F 65(0;)] • 

LJ iO 

In order to prove (3.1), it suffices to prove the following for any conhguration u: 

Buj + i?s(aj) + + Ds(tj) = A{/3Xe) + -As{u)) + . (3-2) 

When 7 ( 0 ;) does not go through any of the diamond edges bordering /, neither does 
7 ( 5 ( 0 ;)). All the contributions then vanish and identity (3.2) trivially holds. Thus we 
may assume that 7 ( 0 ;) passes through at least one edge bordering /. The interface enters 
/ through either A or C and leaves through B or D. Without loss of generality, we 
assume that it enters hrst through A and leaves last through D] the other cases are 
treated similarly. 



Figure 3: Two associated configurations uj and s{uj) 

Two cases can occur (see Figure 3): Either the exploration curve, after arriving 
through A, leaves through B and then returns a second time through C, leaving through 
D] or the exploration curve arrives through A and leaves through D, with B and C 
belonging to a loop. Since the involution exchanges the two cases, we can assume that 
UJ corresponds to the hrst case. Knowing the term A^, it is possible to compute the 
contributions of uj and 5 ( 0 ;) to all of the edges bordering /. Indeed, 

• the probability of s{u) is equal to Pxe-^/q times the probability of u (due to the fact 
that there is one additional loop, and the primal edge crossing / is open); 

• windings of the curve can be expressed using the winding of the edge A. For 
instance, the winding of B in the conhguration u is equal to the winding of the 
edge A plus an additional —6*e turn. 

Contributions are computed in the following table. 


conhguration 

A 

B 

C 

D 

u 


e""=A, 

e-""A, 


5 ( 0 ;) 

I3xe-^A^ 

0 

0 



Using the identity 7 - e ^^2 = we deduce (3.2) by summing the contributions 
of all the edges bordering /. □ 

For a set E of edges of T[G^], deE denotes the set of edges of \ E bordering the 
same face as an edge of E. Also dehne Eint to be the set of edges of the diamond graph 
between two faces of 
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Proposition 11 Fix (3 < 1 and g > 4. Let G^o satisfying (BAPq). There exists Ci = 
Ci{(3,q,9) < oo such that 

He) <c\Y. 

sGE e^deE 


for any E C Eint- 


Proof Sum the identity (3.1) over all faces bordered by an edge in E. It gives a weighted 
sum of E{e) (with coefficients denoted by c(e)) which is identical to zero: 

0 = 5;c(e)f(e)+ 5 ; c(e)F(e). (3,3) 

gG^E e^deE 


For an edge e E E, E{e) will appear in two such identities, corresponding to the two 
faces it borders. Since a loop going through e comes from one of the faces and enters 
through the other one, the coefficients will be —1 and A{f3xe)- Thus, for e E E, E{e) 
enters the sum with a coefficient 


c(e) = A{f3xe) 


\f3xe + e"'^^ 

^ ( 1 ^ (/^ - l)a;. 


1 = 


Xe + e “"a \ / ^Xe + e^'a 
Xe + e”"^ / \/3Xe + 


Xe + e 2 

= W- 


{/3 - l)Xe 
/3xe + e“'^^ 

— e 2 


2(1 — /d)a:eSinh((T|) 


(xe + e'^2)(/ 5 a;e + e ”'2)y (^Xe + 2)(^fjxe + e~ 
> 2(1 — (3) min [xe : e G T[G]} sinh((j|). 


In the second equality, we used that A(xe) = 1. Because of the bounded angle property, 
Xe is bounded away from 0 and oo uniformly, and so is c(e). 

For an edge e G dgE, E{e) will appear in exactly one identity, (corresponding to the 
face that it shares with an edge of E). The coefficient will be A{f3xe) or —1, depending 
on the orientation of e with respect to the face. Thus E{e) will enter the sum with a 
coefficient c(e) which is bounded from above uniformly in e (thanks to the Bounded Angle 
Property). The proposition follows immediately by setting 


max{|c(e)| : e G dgE} max{|c(e)| : e G E[Gl^]} 

min{|c(e)| : e E E} min{|c(e)| : e G E[Gl^]} 


Note that Gi depends only on (3, q and 9, but not on G^o or E. 


□ 


Write dS[G'^] for the edge-boundary of G'^, meaning the set of diamond edges con¬ 
necting two boundary vertices. Note that it is simply E[G^] \ E\^t. 


Lemma 12 Let {G, a, b) be a Dobrushin domain. For a site u on the free arc dab o,nd 
e E dS[G^] a diamond edge incident to u, we have 

F{e) = ^(u EE wired arc dba), 

where IF(e, e^) is the winding of an arbitrary curve on the diamond graph from e to et- 
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Proof Let m be a site of the free arc dab and recall that the exploration path is the 
interface between the open clnster connected to the wired arc and the dnal open clnster 
connected to the free arc. Since u belongs to the free arc, u is connected to the wired arc 
if and only if e belongs to the exploration path, so that 

^ wi’^ed arc dba) = e 7). 

The edge e being on the bonndary, the exploration path cannot wind aronnd it, so that 
the winding of the cnrve is deterministic. Call this winding W{e,eb)- We dednce from 
this remark that 


F(e) = 




tW. 




□ 


We are now in a position to prove onr key proposition. 

Proposition 13 Fix q > A and Goo satisfying (BAPq). There exists C 2 = C2{q,0) < 00 
such that 

^ F(e) < 
e&d£[Go] 

for any (d < 1 and any Dobrushin domain (G, 0,0) (here, 0 is assumed to belong to 
dG). Above, nnd Wmax ore the minimal and maximal winding when going along the 

boundary of G, starting from 0. 

Note that Dobrnshin bonndary conditions on (G, 0,0) coincide with free bonndary 
conditions. 


Proof Let ns start with the case (3 = 1. Snm Identity (3.1) over all . Since A{xe) = 1 
for any e, we find that c(e) = 0 for any e G L^int and c(e) = ±1 on dElG"^] depending 
on the fact that a loop going throngh e points ontwards or inward Bonndary edges 
corresponding to a loop pointing ontward are called exiting, those for which the loop is 
pointing inward entering. We find 

F(e) - = O' 

e exiting e entering 

Since edges exiting or entering belong to dS[G^], Lemma 12 implies that 

F(e) = M 0|, 

for u the site of G bordering e. Note that each vertex u G dG is the end-point of a nniqne 
entering edge, called ein(M) and a nniqne exiting edge eex(n)- With this definition and 
the two previons displayed eqnalities, we find 

[e-^hexW.e,) _ g-7W(ei„(n),e,)j ^^ 0) = 0 

u£dG 

which can be rewritten as 


E f 

«e 9 G\{ 0 } 


,o-Vl/(ee 


_ gO-lV(ei, 


G)Xb) 


T </>G.l,g(^ M 0) = 


^ g<TW(ei. 


(0),eb) _ ^aW{ee 


:(0),e„) 
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Now, when u = 0, ein(O) = and eex(O) = Cf,. Since W{ea, eh) < 27i, 

gCrW(ein(0),e6) _ gCTW(eex(0),ei,) _ ^aW{ea,ei,) _ gO-Vy(ej,,ej,) ^ ^a2n 

When K 7 ^ 0, we find that 


hh(eex(w)) Cfe) hh(ein('u), e^) Vh(eex(w)) Cin(^)) ^ ^[uv] ^ 

where n is a neighbor of u outside of G, and [uv\ is the edge between u and v. Note that 
the existence of this neighbor is guaranteed by the fact that u G dG. Next, 

gCrW(eex(«),ei,) _ QCrW{ein{u),ei,) _ |■g(7W(eex(M),ein(«)) _ ^(7W(ein(u),ei,) y ]^^ae _ j^jgCrWmin 

Therefore, 

[e^e _ i —> 0) < 

«e 9 G\{ 0 } 

Finally, observe that 


E 

F(e)= Y, ^ 0| 


e&d£[Go] 

eed£[G^] 



< 4g<TlVniax Y 0G,1,<7[“ ^^ 0] < 

+ 46""^““ ^ 0G,1,J«^ 


u£dG 

uedG\{o} 


pCT(Wmax~VFrnin) 



^ 1 1 c'^^^ 

- e'^® - 1 



In the first inequality, we used the fact that at most four edges correspond to a boundary 
vertex. 

The case /3 < 1 follows readily since F{e) = Vg/? q[u <—> 0] is an increasing 

quantity in /3. □ 


For a graph G, let us introduce the following, recursively constructed graphs. Let 
(^( 0 ) _ Q Q{k) _ Q{k-i) ^ for any k > 1. They can be seen as successive 

“pealings” of G, each step consisting in removing the boundary of the existing graph. Let 
Ek = Note that Eq = = £[G^] \ d£[G% 

Corollary 14 Let Goo be an isoradial graph satisfying (BAPe). Consider the Dobrushin 
domain (G,0,0) (as above, 0 is assumed to be on the boundary of G). Eor any (d <1 and 
q> A, 

V E{e) < 

eeE, + 



Proof Proposition 11 can be applied to Ek C F^mt to give 


eeSfe eeEkUdeEk 


Since Ek U deEk C Ek-i and F(e) > 0, 


E ^ 

eSEfc 


Cl 

1 + Gi 


E 

eSi?fc_i 
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Using the previous bound iteratively, and Proposition 11 one last time (in the second 
inequality), we hnd 





The claim follows by bounding the sum on the right-hand side using Proposition 13. □ 


The study above can be performed with F instead of F. We obtain the following 
corollary. 


Corollary 15 Let Goo satisfying (BAPg). Consider the Dobrushin domain (G, 0, 0) (0 
is assumed to be on the boundary of G). For any fd < 1 and q> 


V F(e) < 

eeE, + 


4 Proof of Theorem 2 

Without loss of generality, we assume that n = 0. Fix/3<1. We aim to prove that there 
exists c = c{/3, q,9) > 0 such that 

^^ ^ exp(-c|M|) 

for any u E Goo containing 0. We now hx Goo containing 0 and satisfying (BAP^). We 
stress out that the constants involved in the proof depend on 9 only but not on Goo or u. 

The case g = 4 is derived through stochastic domination between random cluster 
measures. Indeed, for every /3 < 1, there exists (/?', q) with g > 4 and /3' < 1 such that 
the random-cluster measure 4>g^ y q stochastically dominates the random-cluster measure 
*?^Goo/ 3 4' refer to [27, Theorem (3.23)] for details on this fact. It follows from this 
stochastic domination that 

‘/’Goo,/ 3 , 4(0 ^^ ^^ “) < exp [ - c(/?', g, 9) |m|] 

for any u E Goo- It is therefore sufficient to assume that g > 4, which we now do. 

For a graph Goo, we identify a convex subset A of with the subgraph given by 
vertices in Goo H A and edges between them. The set dA is referring to d{Goo H A). Note 
that this set is not necessarily connected, but this will not be relevant in the following. 
For r > 0, let i?o(r) = {a: G : |a;| < r}. 

Lemma 16 There exists ci = ci(/3, g, 0) >0. Assume that v E Goo is on the positive real 
axis. Then, 

</>[0,oo)x(-oo,oo),/3,g(0 ^^ '^) ^ exp[-ci|n|]. 

There could be no vertex v on the positive axis, but there is no loss of generality in 
assuming that v belongs to this positive axis, since the graph Goo can be rotated around 
the origin in order to obtain an estimate valid for any vertex v E Goo, where the half-plane 
[0, 00 ) X (— 00 , 00 ) is replaced by a half-plane containing v whose boundary contains 0 
and is orthogonal to the vector given by the coordinates of v. 
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Proof Define Gm,n to be the connected component of the origin in [0,n] x [—m,m]. Set 
k = \v\/2. Let m > 100 and n > max{2|n|, 100} (these two conditions avoid nasty local 
problems). Apply Corollaries 14 and 15 to Gm,n and k to hnd (we use the notation of 
the corollary) 


F(e) + F(e) < 


Cl 


1 + Ci 


k 


The maximum and minimum windings on dE[Gm,ri\ are bounded by a certain constant 
Cs < oo. This statement comes from the fact that the winding of dE[Gm,'n\ is roughly 
comparable of the winding of the boundary of [0,n] x [—m, m] (there are a few local 
effects to take care of). Let us sketch an argument. There exists a path of adjacent faces 
of the subgraph = [0, n] x [—m, m] \ [4, n — 4] x [—m + 4, m — 4] of Gm,n going 
around [4, n — 4] x [—m + 4, m — 4]. The constant 4 has been chosen to £t our purpose. 
It can probably be improved but this would be of no interest for the proof. In particular, 
dE[Gm,n\ is contained in Hm,n- The bounded angle property shows that two vertices 
cannot be arbitrary close to each other, which prevents the existence of paths of edges 
in E[Hm,n] winding arbitrary often around a point of As a consequence, the winding 
along dE[Gm^n] is indeed bounded by a universal constant. 

The previous bounds on ITmax and ITmin imply 


^F{e) + F(e)<cJ-^ 

.EE. + 


Next, observe that 

F{e) + F{e) = e 7 ). 


Let Co be the cluster of the origin and c^outCo its outer boundary, meaning the set of 
vertices connected to the free arc by a path in E[Goo] \ E[Co\. A diamond edge e = [vy], 
where v G Gm,n and y G G*^^ belongs to 7 if and only if 0 is connected to v and y is 
connected to the free arc. In other words, v is on the outer boundary of the cluster of 0. 
We deduce that 


•#-L,.A,(9o„.c„nGW T-ejs 5^ 6 sc„) 

m,n 

< 5 ^ FW + 7e) 


eeE> 


(fc-l) 


< C 4 


C'l 


I + Cl 


fc-1 


Letting m go to inhnity and using the uniform bound above, 

n gW„ 7 ^ 0) < c, j , 

where Goo,n = [0,n] x (— 00 , cx)). Next, observe that Goo.n does not contain any infinite 
cluster 0 G'^„,/ 3 ,q-almost surely (in fact, this is true for any (3' < 00 , since the graph is 
rough isometric to Z). Let us provide a rigorous proof of this statement. The experienced 
reader can skip the next paragraph. 
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The bounded angle condition implies that the weight x^, is bounded from above uni¬ 
formly on Goo- Therefore, there exits C 2 = C 2 (/ 3 , q,9) > {] such that for any hnite set S of 
dual edges of cardinality r 


'^Goo.n,/ 3 ,<?(^very edge in S is dual-open | > exp[-C 2 r]. 


The constant C 2 comes from the hnite-energy of the random-cluster model, i. e. the prop¬ 
erty that the probability for an edge to be closed is bounded away from 0 uniformly in 
the state of all the other edges; see [27, Equation (3.4)]. Next, it is possible to divide 
the strip into an inhnite number of hnite pieces by considering disconnecting paths P of 
length less than c^n for some constant C 3 = 03 ( 6 ^) < 00 . The existence of these paths 
is easily proved by considering a sequence of set of adjacent faces cutting the strip, and 
by noticing that the bounded angle property bounds from above the number of edges 
bordering a face by vr/sin(6'/2). Now, conditioned on the state of the others paths, each 
one of these paths has probability larger than exp[—C 2 C 3 n] > 0 of being dual-open. This 
immediately implies that there is no inhnite cluster (pGaa surely. 

Since there is no inhnite cluster almost surely, Cq intersects G^]n if and only if c^outCo 
intersects Furthermore, edges have length smaller than 2, which implies that 

U e Gt]n (we use the fact that n — |u| and |u| are larger than 2k). Hence, 


^ < 0L,.,/3,,(ConGg„ ^ 0) = 0^^ ,^,^,,(5ConGg„ ^ 0) < Gs 


C'l 


1 + Gi 


The proof follows by letting n go to inhnity and then by choosing Ci = Ci{f3,q,6) > 0 
small enough. □ 



Figure 4: The gray area is 7 . The dual path surrounding it is T. It is the 
exterior-most dual circuit. 

We are now in a position to prove Theorem 2. Let n > |u| -|- 2. We work with 
the random-cluster measure on Bo{n) with free boundary conditions. Let X^ax be the 
site of Bo{n) n Cq which maximizes its Euclidean distance to the origin (when several 
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such sites exist, take the first one for an arbitrary indexation of sites in Goo)- Note that 
\u\ < |-^max| < n. Therefore, 

^^ «) < ^ \ Boi\u\) : = v) 

— 4 ’%o{n),l 3 ,qi^niax = v). ( 4 . 1 ) 

For V, let (t(n) be the set of dual-open self-avoiding circuits 7 surrounding the origin and 
V and such that any site of surrounded by 7 is in Bq{\v\). Let 7 be the set of sites 

of Boin) surrounded by 7 G C(n), see Fig. 4 . Dual-open circuits in <t{y) are naturally 
ordered via the following order relation; 7 is more exterior than 7' if 7 C 7'. 

If Xmax = 'y, then V is connected to 0 and there exists a circuit in £(n) which is 
dual-open (simply take the boundary of Cq, drawn on the dual graph). Let {F = 7} 
be the event that 7 is the exterior-most dual-open circuit in £(n). With this definition, 
X ma v = n if and only if v is connected to 0 and there exists 7 G (t(n) such that F = 7. 
Therefore 


0 Bo(n),/ 3 ,g(^max = v) = (/)Bo(n),/^,q (0 < —> V and 37 G ^(n) i F = 7) 

= (/>Bo(n)Ag (0 M n and F = 7) 

7e£(i)) 

= ^ ^|r = 7)</>Bo(n),/3,,(r = l) 

'yG£(v) 

Since F is the exterior most circuit in £(n), {F = 7} is measurable with respect to 
dual-edges outside 7. Furthermore, edges of 7 are dual-open on {F = 7}, see Fig. 4. 
Therefore, conditioned on {F = 7}, the measure inside 7 is a random-cluster model with 
free boundary conditions. Hence, 

<^Bo(n)Ag(0 ^ = 7 ) = 07 / 3 ,g (0 ^1 u). 

Rotate the graph Bo{n) in such a way that v is on the positive axis. Let H = {— 00 , n] x M. 
Observe that by definition of £(n), ^ G H. The comparison between boundary conditions 
leads to 


07 / 3 , q (0 ^— >v) < 0?^,/3,g(O ^t n in 7) < i — > v) < exp[-ci|n|] 

by Lemma 16 . This implies 

0BoW,/3,</(^max = t') < exp[-ci,j(F = 7) < exp[-ci|n|]. 

7ec(//) 

In the second equality, we used the fact that the union of {F = 7} for 7 G £(n) is disjoint. 
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Going back to (4.1), we find 


exp[-ci|n|] 

'ueBo(n)\Bo(|M|) 

Y exp[-ci|n|] 

■ueGoolSodul) 

< Y^ \Bo{\u\ + k + 1)\ Bo{\u\ + k)\ exp[-ci(|M| + k)] 

k>0 

^ + ^ + 1)1 exp [-Cl (|m I + k)] 

k>0 

< C 4 |m| exp[— ci|m|]. 

In the fourth inequality, we used the fact that the number of sites in i?o(c + 1) \ 5o(c) is 
bounded by sip(g/ 2 ) 1) because of the bounded angle property. 

The proof follows by letting n go to infinity and by choosing c = c(/9, g, 0) > 0 small 
enough. 

Remark 17 The main point of the previous proofs, and of our whole approach, is that 
the observable F behaves like a massive-harmonic function. ITe are not able to make 
exact computations though, because we made no particular effort to choose our coupling 
constants with this goal in mind (all that matters here is that the constant Ci is bounded). 
In a recent work [9], Boutillier, de Tiliere and Raschel define a massive Laplacian operator 
on isoradial graphs for which they get the explicit rate of exponential decay of massive- 
harmonic functions. It would be interesting to see whether a better tuned choice of the 
Peifi) can lead to an exactly massive-harmonic observable in their sense, as this would 
lead to sharper asymptotics for the two-point function of the model in that case. 

5 Proofs of the remaining theorems and corollaries 

Proof of Theorem 4 Fix /3 < 1. Let Goo be an infinite isoradial graph. Without 
loss of generality (simply translate the graph), we assume that ?/ = 0 G G^o- For r > 0, 
Theorem 2 implies 

0L,/3,g(O ^ dBo{r)) < Y ^ “) 

u&dBo{r) 

7rr‘^ 

-4;R^“p[-c(/3.9.«)(’'-2)]- (5-1) 

In the second inequality, we used the fact that any site u G dBfir) is at distance larger 
than r — 2 of the origin since any primal edge is of length smaller than 2 (the diameter 
of circles is 2), and that u G OBfir) is connected to a site outside Ro(^)- We also used 
the fact that the cardinality of i?o(^) is smaller than 7rr^/(4sin(0/2)) since any edge of 
Goo corresponds to a face of of volume larger than 4sin(6'/2) thanks to the bounded 
angle property. Letting r go to infinity, we obtain that ficao,i3,q^^ ^ ~ b- 

Let us now consider with (3 > 1. For a dual vertex y E G*, let A{y) be the event 
that there exists a dual-open circuit surrounding the origin, i.e. a path of dual-open edges 
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disconnecting 0 from infinity in Observe that the dual circuit must go to distance 
\y\. Since the dual model is a random-cluster model on the dual isoradial graph, with 
free boundary conditions and with f3* = 1/f3 < 1, (5.1) implies 

2 

for any y G G*^. Borel-Cantelli lemma together with the bound O < 4 gjp(g/ 2 ) 

implies that 

</>Goo,/ 3 ,g(^here exist inhnitely many y e G*^ : A{y)) = 0. 

This immediately implies that there exists an inhnite cluster (pQ^ ^ ^-almost surely. □ 


Proof of Theorem 3 Let us now turn to the uniqueness question. This is the only 
place where a priori information on uniqueness is used. Recall that in this 
case, T)c^,q dehned in Proposition 8 is countable. Since (pGao,y,q ~ ^Goa,i3,g outside of the 
countable set T>G^,q, for any (3 < 1 there exists (3 < (3' < 1 such that 4>Goo,P’,q ~ 

We deduce that 

M OO) < 0L,/3',g(O ^ 1 OO) = 0L,;3/,g(O ^^ CX)) = 0. 

From [27, Theorem (5.33)], we deduce that 0 Gcx>,/3,<? ~ ^G^,i3,q any (3 < 1. Theorem 
(5.33) is proved in the case of but the proof extends to the context of periodic graphs 
mutatis mutandis. By duality, the inhnite-volume measure is unique except possibly for 
(3 = 1. This implies in particular that there is ^ ^-almost surely an inhnite cluster 
whenever (3 > 1. □ 


Proof of Corollary 5 We present the proof in the case of the square lattice, the cases 
of triangular and hexagonal lattices being the same. If then XiX 2 = 1 

— - By embedding the square lattice in such a way that every face 


where Xi = 


is a rectangle inscribed in a circle of radius 1, and the aspect ratio is given by X\lx 2 , 
we obtain an isoradial graph with critical weights Xi and X 2 . The model is therefore 
subcritical (respectively supercritical) if and only il (3 <1 (respectively (3 > 1). □ 


Proof of Corollary 6 Corollary 6 follows directly from the classical coupling between 
random-cluster models and Potts models. □ 
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